Long- wavelength approximation for string 
cosmology with barotropic perfect fluid 



Piret Kuusk 1 

Institute of Physics, University of Tartu, Riia 142, Tartu 51014, Estonia 

Margus Saal 2 
University of Tartu, Tahe 4, Tartu 51010, Estonia 

PACS number: 98.80 Dr 



Abstract 

The field equations derived from the low energy string effective ac- 
tion with a matter tensor describing a perfect fluid with a barotropic 
equation of state are solved iteratively using the long-wavelength ap- 
proximation, i.e. the field equations are expanded by the number of 
spatial gradients. In the zero order, a quasi-isotropic solution is pre- 
sented and compared with the general solution of the pure dilaton 
gravity. Possible cosmological models are analyzed from the point of 
view of the pre-big bang scenario. The second order solutions are found 
and their growing and decaying parts are studied. 
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1 Introduction 



The long-wavelength iteration scheme for investigating the early stages of 
cosmological models was first introduced by Lifschitz and Khalatnikov [1] in 
the case of the Einstein equations with pure radiation as a source term. An 
improved scheme was presented by Tomita [2] where its main assumption 
was taken to be the smallness of spatial variations of metric dijij compared 
with time variations <9o7rr Then Comer, Deruelle, Langlois, and Parry [3] 
applied it to the Einstein equations with a perfect fluid matter (with a 
barotropic equation of state, cf also [4]) and in the case of a scalar field 
with a potential. The next step was done by Tomita and Deruelle [5] who 
examined the case of two fluids (an inflationary and a standard fluid). Fi- 
nally, the scalar-tensor theory with a perfect fluid has been discussed by 
Comer, Deruelle, and Langlois [6]. The long- wavelength approximation for 
the Hamilton-Jacobi formalism (the gradient expansion for the generating 
functional) was developed by Salopek, Stewart, and Parry [7]. 

In the following, we solve iteratively field equations derived from the low 
energy string effective action with a perfect fluid matter with additional re- 
quirement that the matter obeys a barotropic equation of state. In general, 
the approximation scheme can be applied in the case of anisotropic and inho- 
mogeneous spacetimes. In the present paper we ignore the local anisotropy, 
although this assumption may be incorrect near singularity [8]. We find a 
quasi-isotropic solution for the field equations in the zero approximation (ne- 
glecting all spatial gradients) which include both post-big bang and pre-big 
bang [9, 10, 11, 12] cosmological models. It has been demonstrated [13, 14] 
that a smooth transition from a pre-big bang branch to a post-big bang 
branch (the graceful exit problem, see also [15, 16, 17]) is impossible in this 
approximation, even if we incorporate some dilaton potential, or cosmolog- 
ical constant. We analyze the nature of discontinuities which are contained 
in cosmological models based on a quasi-isotropic solution with a barotropic 
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perfect fluid matter source (cf [18]). 

Recently, Khoury, Ovrut, Seiberg, Steinhardt, and Turok [19] presented 
a new cosmological scenario where the Universe is initially contracting to- 
wards a big crunch and then makes a transition through a singularity to 
the post-big bang Universe. The model is a part of the M-theory motivated 
ekpyrotic scenario developed by Khoury, Ovrut, Steinhardt, and Turok [20] 
and occurs in a recent proposal which introduces the cyclic nature of the 
ekpyrotic Universe [21]. The same considerations can be applied to the 
reversal problem in the pre-big bang scenario. 

The paper is organized as follows. In the next section, the string frame 
field equations in a synchronous coordinate system are given. In the third 
section, solutions of the field equations in the zero approximation (taking 
spatial gradient terms to vanish) in the case of pure dilaton gravity and in 
the case of dilaton gravity with a perfect fluid matter obeying a barotropic 
equation of state are presented. Specific interpretations of solutions are 
considered and problems of cosmological model building are studied. In the 
fourth section, the second order solutions in the long-wavelength approxima- 
tion are found. It is analyzed how the equation of state for matter influences 
the evolution of inhomogeneities. The fifth section is a summary. 

2 Field Equations 

Upon compactification the low energy action can be written in the lowest 
order in the inverse string tension and coupling and using the string (Jordan) 
frame as 



<t>\ 4. 



K + g^d^dA - ±-H^ p H^p + V{<f>)\ 




(1) 
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Here 4 TZ is the curvature scalar of the metric g^ u , (f) is the dilaton field 
determining the strength of the gravitational coupling through g 2 ~ e^, 
H^ up is an antisymmetric field strength, X s is the fundamental string length 
scale, L m is the matter lagrangian, and A is the cosmological constant. We 
assume that there is no direct coupling between the dilaton and matter 
fields. 

The energy-momentum tensor for a perfect fluid matter is 

Tp V = + p)u^u v + pg^ (2) 

and we assume that the energy density and pressure are related by a 
barotropic equation of state 

p = (r - i) P , o < r < 2. (3) 

Here r = | corresponds to a pure radiation, T = 1 to a pressureless dust, 
r = | to an inflationary fluid, T = to a phenomenological matter sector 
cosmological constant, and T = 2 to a stiff fluid. 
Let us introduce a synchronous gauge 

ds 2 = -dt 2 + jij (t, x)dx i dx j (4) 

and additional simplifying assumptions of vanishing three- form field, H^ up = 
0, dilaton potential, V(<j>) = 0, and cosmological constant, A = 0. Now the 
field equations following from the action functional (1) with a barotropic 
perfect fluid (2), (3) as the matter source can be written 

nl = \k+-K)K) = \^> P {2 - 3r - 2Tu s u s ) + 

I 4 J J 4 4 

+ ^(30-0 2 + cU<9 r 0- V r V», (5) 
K^hpiK-VjKi) = ieVWl + u s u° - 1 -K{v ^ + ,(6) 
ni = R{ + ±ki + \KKi = le*p[{2-T)8i + 2Tu iU i]+ l -Kii> 
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+^Kj4 + U{{4> -<p 2 + d r 4>d r 4> - v r v» - ViV^ , 



(7) 



4> + ^K4> -ft- V r V> + d r 4>d r 4> = ^e*[(3r - 4)p] . (8) 

Here = 7^, K = K\ = j^jij and is the Ricci tensor of the three- 
metric 7y. In addition, we have the usual conservation law for the matter 
density V M T^ = 0. 

3 Zero Order Solutions 

In the zero order approximation we ignore spatial gradients and local 
anisotropy in field equations (5)— (8). A solution for the metric is taken 
to be in a quasi-isotropic form 

V j (t,x k )=a 2 (t)h ij (x k ), (9) 

where h{j(x k ) is a time independent seed metric and a(t) is a scale fac- 
tor of an isotropic and homogeneous Universe. As demonstrated in [2, 3], 
the quasi-isotropic solution (9) can be regarded as an attractor of general 
anisotropic solutions. We will consider only a particular class of solutions 
with a power law evolution of the scale factor and a logarithmically evolving 
dilaton 

a = a Q T a{T \ = o -/?(r)lnr, (10) 

where r is a time parameter. Solutions of this form can be found only in the 
case of the pure dilaton gravity with a vanishing potential, V{4>) = 0, and 
with a perfect barotropic fluid matter (2), (3) as a source. Inclusion of the 
cosmological constant A introduces a hyperbolic nature of solutions and we 
shall not consider this case here. Solutions (10) inevitably contain singular 
points (a — > 0, g 2 ~ — ► 00) and it is not possible to avoid them. 
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3.1 Pure Dilaton Gravity 



The field equations for the Hubble parameter H(t) = 
4>(t) can be written as follows 



and the dilaton 



H + 3H 2 -H<j> = 0, 
6H 2 - 6H<j) + <j) 2 = 0, 
4> + 3H<p - (f 2 = 0. 



(11) 
(12) 

(13) 



The system is not overdetermined, since only two equations are independent. 
The most general isotropic solution reads 



Here B = ±1 denote two branches of solutions. The range of variation of the 
time coordinate t must be chosen such that At + to > 0, so the time inversion 
amounts to a change of constant A — ► —A. The dilaton gravity analogue of 
a post-big bang model corresponds to A > 0, B = 1, a standard pre-big bang 
model [9, 10, 11, 12] corresponds to A < 0, B = —1. Integration constants 
ao, 4>o are equal to values of a(t), 4>(t) at an initial moment U = and A 
determines initial values for a(t), 4>{t). 

An one-to-one correspondence between post-big bang solutions and pre- 
big bang solutions can be arranged by duality transformations 



a(t) -► o(t) = a^i-t), (f)(t)^(t){t) = <f)(-t)-6\na(-t). (16) 



In terms of general solution (14), (15), the correspondence between a partic- 
ular post-big bang solution (B = 1) with ao, 4>o and a pre-big bang solution 
(B = —1) is encoded in integration constants ao,^o,A that have undergone 
the same duality transformation (16): 



B 



a = a (At + t )Vs } 
(j) = ci)o + (BV3-l)ln(At + to). 



(14) 



(15) 



a(t) = ao 1 (t - At) 
$(t) = (j) - 6 In a - (VS + 1) ln(i - At). 



(17) 



(18) 



6 



3.2 Dilaton Gravity and Perfect Fluid Matter 

After some manipulations the field equations become 

H + 3H 2 -Hj> = ^p(r-l), (19) 

6H 2 - 6Hj> + <j) 2 = e^p , (20) 

+ 3^ _ ft = I e ^( 3 r - 4) , (21) 

p + 3i3T> = , (22) 

where the last equation is the conservation law of matter. Loosely speaking, 

equations (19), (21) and (22) describe the evolution of H, cf> and p, remaining 
equation (20) imposes a constraint. 

A quasi-isotropic solution of the field equations can be written as follows 

a = a (t ±1x1*)^, (23) 

<t> = 0o~ 2(4 A 3r) M*o±lxl*), (24) 
6r(r-i) 

p = p (t ±\x\t) — , (25) 

where 



_ e*°po (A) 2 
X = ± i -4 (A-2) ' (26) 



A = 3(r - l) 2 + 1 > 0, (27) 

and ao, <Pq, po are integration constants which are equal to the values of 
functions a(t) , <p(t) , p(t) at an initial moment U = ^j^- The solution is 
not the most general one, since initial values of a(t) and <p(t) are both 
uniquely determined by barotropic index T and the initial values of <fi(t) 
and p(t); in the most general case they must contain an additional arbitrary 
integration constant. As a result, solution (23)-(27) turns out to be singular 
at T s = 1 ± where A — 2 = 0, x —> 00 • Note that in the case of 
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unusal matter with (T — l) 2 > 1/3 we must take po < 0. For the zero 
order three-velocity Uj, field equation (6) implies Ui = 0. The range of 
variation of the time coordinate t must be chosen such that to ± Ixl* ^ 0; 
to is a numerical constant corresponding to the freedom of constant shift in 
time. In comparison with solutions of the pure dilaton gravity (14)-(15), 
the integration constant A is replaced by constant x = xiPOi 'Po-, F) and the 
choice of the branch, B = ±1, is imitated by the choice of barotropic index 
T, or more exactly, by the sign of expression T — 1 € [—1, +1]. 
Solution (23) for a(t) determines the Hubble parameter as 

±IyI 2 • -y 2 2 

h= , m — 4(r- 1), h = - — ^L_4(r-i). (28) 

(to ± Ixl*) A ^ ; ' (t ±\ x \t) 2 A y ' y ' 

We see that constant x is proportional to the initial value of the Hubble 
parameter at the initial moment tj, to ± \x\U = lj an d in the case of matter 
with barotropic index r s = 1 ± -i it diverges. In what follows we don't 
consider barotropic indices belonging to the neighborhood of T s . 

If we take T > 1 and temporal argument (to + \x\t) we get a decelerating 
post-big bang Universe, e.g. for T = | (radiation dominated stage) we 
have \x\ = \J\ e ^°Po an d a (t) = a o(*o + Ixl*) 5 ) matter density is decreasing, 
p = po(to + |x|i)~ 2 , and dilaton is freezed to a constant, (p = 4>o- Integration 
constants ao and po are (initial) values of corresponding variables at the 
moment tj = ^j 1 - A singularity (a — > 0, p — > oo) is reached at the moment 
t s = — 1^| • The domain of the time variable is t G (t s ,+oo). Taking into 
account our choice of initial moment tj we see that U = t s + j^-, i.e. ^ is 
the time interval from the singularity to the initial moment. If to > we 
have t s < U < and the Universe is regular at the moment t = 0. If the 
source matter is exotic (r < 1), a solution (23)-(25) represents a contracting 
Universe beginning from a singularity at t s and going through tj. 

The pre-big bang branch of solutions corresponds to the minus sign in 
front of \x\t and exotic matter T — 1 < 0. The domain of the time coordinate 
is t G (— oo, = too), i-e. to — \x\t > 0, and the behaviour of the scale factor 
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is super inflationary (d > 0, a > 0, H > 0, H > 0). Integration constants ao, 
00, and po are the values of corresponding variables at the moment ij = 
and here ij = t s — i.e. t^t is the time interval from the initial moment to 
the singularity. A contracting solution with a future singularity is obtained 
if r > 1. 

Expressions (24) for (j) and (28) for H imply 

H = -(r-i)4>, 4> = (j)-w (29) 

and the behaviour of solutions depending on barotropic index V can be 
summarized on a phase diagram analogous to the corresponding one familiar 
from the pre-big bang scenario of the pure dilaton gravity: 




r- 1 <o 
r-i >o 



to + btlt =(0,°o) t -lxl t =(0,oo) 

te(t S) ao) tec-oo.t^) 



As in the case of pure dilaton gravity, solution (23)-(25) contains both 
the post-big bang and the pre-big bang branch. The latter ones can be 
related to the former ones by duality transformations (16), but in order to 
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satisfy the field equations, here we need an additional transformation for 
barotropic index T and for energy density p: 

= =>f->2-r, p^pa 6 (-t). (30) 

P P P 

The corresponding pre-big bang solutions (i.e. with the time dependence in 
the form (to — |x|i)) now read 

a = ao 1 (t - IxK) 2 ^ , (31) 
= <ft -61na - 2 ( 4 ~ 3r ) ln(t -| x |f), (32) 

„ 6f(f-l) 

P = Poag (*o - Ixl*) ~- (33) 

Note that poe^° = poe^°, A = A and to don't change under duality trans- 
formations (16), (30). 

3.3 Model building 

In the case of pure dilaton gravity, the main idea of the original pre-big 
bang cosmology was to use a superinflationary pre-big bang branch obtained 
from a decelerating post-big bang branch by duality transformations for 
constructing a cosmological model valid in the time range t £ (— oo,+oo). 
However, due to rather general no-go theorems [14] a smooth connection of 
these two branches turned out to be impossible, at least without some kind 
of higher order corrections. In the case of a perfect barotropic fluid, mutually 
dual solutions (23)-(25) and (31)-(33) contain different barotropic indices 
T, f = 2 — T and consequently also different matter sources. This means 
that if we try to repeat the strategy of the original pre-big bang cosmology, 
we must admit a profound change in the type of cosmological matter and 
join a pre-big bang solution with barotropic index (2 — T) with a solution 
with barotropic index T (cf [17, 18]). We are free to choose the moment of 
branch change, but the explicit solution (28) for H indicates that we cannot 
achieve continuity for both H and H. 
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For example, let ao = 1 and to = 1 and let the junction moment be the 
regular point t = 0. From the solutions (23)— (25), (31)-(33) it is easy to see 
that a(t), a(t), (f>{t) and p(t) are continuous there, but not H(t) (or a(t)), 
<j)(t) and p(t): 

B(0) = 2 X 2 (T - 1)A- X = -H(0), 0(0) = 0(0) - 4IXIA- 1 , (34) 

M0)=p(0) + 12p |xl(r-l) 2 A- 1 . (35) 

It is possible, that these kinematical discontinuities could be given in terms 
of physical changes of matter sources. However, a dynamical description of 
these changes remains open. Numerical values of discontinuities depend on 
the numerical value of the constant x given in terms of barotropic index T 
and initial values of the dilaton and matter density (26) . We can choose ^ ~ 
-^Planck! then discontinuities are nearly vanishing, but initial moments are 
almost in the singularities and higher order corrections must be taken into 
account. Alternatively, we can choose initial moments to be in a classically 
regular region, far from the singularities; then, however, discontinuities are 
finite. 

The other possibility is to abandon duality transformations altogether 
and to connect solutions with the same barotropic index T, but with different 
temporal arguments (to ± \x\t) an d with arbitrary integration constants ao, 
0o> Po, X an d to (cf ekpyrotic models [19, 20]). However, according to the 
solution (28) if H is continuous, H is not. 

As already mentioned, the moment of branch change can be chosen ar- 
bitrarily and it may be a regular or a singular point of the solution. This 
means that the problem of graceful exit is distinct from the problem of sin- 
gularity. The graceful exit problem in the string cosmology is in fact very 
similar to that of ordinary inflationary cosmology, where accelerating infla- 
tionary stage must be joined with a decelerating FRW model. In ekpyrotic 
models, discontinuities are concealed in the singularity. 
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4 Second Order Solutions 



Let us we adopt the procedure presented in [3, 6] to equations (5)-(8) derived 
from the string effective action (1). The second order solutions for the three- 
metric, dilaton, energy density and three-velocity are assumed to be in the 
form 



Here, Rij(h k i) and R(h k i) are the 3-dimensional Ricci tensor and the scalar 
curvature formed from the seed metric h k i, respectively. The zero order 
solutions a(t), 4>(t), and p(t) are given in the previous section. Substituting 
them into field equations, keeping only terms up to the second order in 
spatial gradients, collecting coefficients at Rj, and R we get a system of 
ordinary differential equations for the second order corrections f 2 , g 2 , and 

<t>2- 

4.1 Pure Dilaton Gravity 

If we substitute expressions (36)-(37) into field equations (5)— (8) and take 
p = we get equations for the second order corrections f2(t), g2(t), and 



jijit, x k ) = a 2 {t) hij + f 2 (t)R tj (h kl ) + -(g 2 (t) - f 2 (t))R(h kl )h ij , (36) 



<P(t,x k ) 
p(t,x k ) 

Ui(t,X k ) 



<f>{t)+<h(t)R(h k i) : 
p(t) + p 2 (t)R(h kl ), 
u(t)ViR(h kl ). 



(37) 
(38) 
(39) 



h + (3ff - 4>)h = ~2a 



-2 



(40) 



(2H - <f>)g 2 - {6H - 2<P)<p 2 = -a 



-2 



(41) 




(42) 
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Upon substituting zero order solutions (23)— (25) into equations (40)— (42) 
we get the second order solutions (we have taken A = 1) 

/±(i) = c f + 4 info ± f ) - ^ * || ± f ) 2 ^ ; (43) 



^ ef + ^4(^r 1 + i ^^(^^, (44) 

^(*) = ^+4fa±r 1 - 2 ^^ 7) fa±^- (45) 

The first terms in expressions (43)-(45) are solutions of corresponding ho- 
mogeneous equations and the last ones are particular solutions of inhomo- 
geneous equations. 

As indicated by Comer, Deruelle and Langlois [6] terms proportional to 
(to it)" 1 which diverge at (to±t) — > can be removed by an infinitesimal co- 
ordinate transformation that preserves the synchronous coordinate system. 
In the second approximation, the transformation of the three-metric can be 
given as 7^ = 7^ — 2aaT2R(hij)hij and the transformation of the dilaton 
as (/)' = (f) — (PT2R with T2 = const [6]. If we take T2 = T^m^ then terms 
proportional to d^(to ± i) -1 cancel. The logarithmic term in fa cannot be 
removed by a coordinate transformation (an analogous problem occurs in 
the case of linear perturbation theory [22]), but probable it can be treated 
by means of a renormalization procedure [23]. 

In the case of lower signs (pre-big bang) solution (43)— (45) coincides with 
the solution of the corresponding Hamilton-Jacobi equation in the second 
approximation presented by us earlier and used for estimating the size of an 
initial homogeneous domain for getting enough inflation in the pre-big bang 
stage [24]. 
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4.2 Dilaton gravity and perfect fluid matter 

Following the procedure described above and eliminating the matter density 
p from field equations we get a system of equations for f2(t), 92{t), and fait) 

h + (3ff - fah = -2a" 2 , (46) 



h + 



6H- 



(4 - 3r) v 



92 ~ 



6(7 -6r)„ 12(1 -r). 



h + 



2H - 



(4 - 3r) 



4-3r 

4-3T r 
6 



92 ~ 



H - 



(4 - 3r) 

-2a" 2 , 
4 



(47) 



4-3r" (4-3r) 
(io -er), =Q 



(48) 



4-3r 

These equations are apparently singular at T = | but in fact they are not 
because of the form of the solution in the zero approximation. 

Let us substitute the first order solutions (23)-(25) into equations (46)- 
(48). Upon long but straightforward calculations we obtain the following 
system of equations 

2 



ff ± ^\x\(to±\x\t)- 1 f? = -2a^it ±\ X \ty 



4(r-i) 



£ ± 2\x\ito±\x\t)-'g 2 t = {ZT-b)av Z it Q ±\ X \t) s 
2 4 — 3r 



4(r-i) 



(49) 



(50) 



(51) 



(5 - 3r)A 

The upper sign corresponds to the post-big bang behaviour and the lower 
sign corresponds to the pre-big bang behaviour. Retaining all integration 
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constants general solutions of equations (49)-(51) can be written 



ft(t) = cf + 4(t ±\x\t)^ 



g£(t) 



3r z -6T+2 



A 2 ixrv 



4(r-i) 



(3r 2 - sr + 6)(3r 2 - lor + 10) 
df + 4(t ±\ x \t)- 1 

(6r-10)A 2 | X |- 2 a 2 

4(3r 2 - sr + 6)(9r 2 - 22r + 16) 

i i , , . . , . 3r 2 -er+2 3r — 4 



(t ±\x\ty-— 



(52) 



6(r - 1 



(t ±\x\t) 2 ~ iJ ^ R , (53) 



■d±(t ±\ X \t)- 1 



2„-2 



+ 



(3r - 4)(3r 2 - 6r + 6)A 2 \ X \~ 2 a 



4(3r 2 - sr + 6)(9r 2 - 22r + i6)(3r 2 - lor + io) 

For the matter density correction pi we get from equation (8) 
and for the three-velocity we get from equation (6) 



(t ±|x|t) 2 -^.(54) 



2 + 3H<p 2 + ^g 2 <j) + 2# 2 - i(3r - 4)e 



u(t) 



2e" 



7^2 - -7^/2 - #</>2 - </>2 



For the density contrast we can write the expression as follows 



S = £ 



p (3r - 4)poe^ 



(t ± \ x \t) 2 x 



6 2 + 3#0 2 + ^ + 2<^> 2 - 2 ( 3r - 4)e*p0 2 



(55) 



(56) 



(57) 



Let us consider the behaviour of inhomogeneities in the pre-big bang 
case near the final singularity (\x\t — ► to)- Terms proportional to (to — 

|x|t) 3 are decaying if T < 1 - and T > 1 + e.g. at T = 

? 4(r-i) 

and T = 2. Terms proportional (to — |xK) A are decaying at all values 
of r. Terms proportional to (to — |xl*) _1 can be treated as in the case 
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of pure dilaton gravity: they can be cancelled by an analogous coordinate 
transformation with T2 = ± i2(r-i)M • 

In the post-big bang case, the model with r = | should describe a radi- 
ation dominated FRW cosmology and so it does in the zero approximation. 
However, in the second approximation the initial inhomogeneities are not 
decaying at late times f > to, although logarithmically divergent terms that 
occur in the solution for the pure dilaton gravity are absent here. The full ex- 
pression for the metric (36) includes terms proportional to (to + \ x\t) 2 which 
are growing at all values of T. An analogous term also appears in solutions 
of field equations corresponding to other models [3], [6], [4]. Khalatnikov, 
Kamenshchik, and Starobinsky [4] have recently indicated its inevitability 
near the singularity. However, in our solution and in the solutions presented 
in [3] it also appears in the late time evolution. 

5 Summary 

In this paper we investigated cosmological solutions of equations (5)-(8), 
derived from the low energy effective string action, in the framework of 
long-wavelength approximation. The zero order solution (23)-(27) contains 
both the pre-big bang and the post-big bang branch, but as it is well known, 
it is not possible to connect these two branches smoothly. We have found, 
that the problem of singularity and the problem of graceful exit may be 
disconnected. The graceful exit problem resembles the problem already 
known from the first models of inflation: how to get a smooth transition from 
an accelerating background to a decelerating background. The transition 
moment and the moment when classical solutions become singular may be 
different, in general. On the other hand, as indicated in the context of 
ekpyrotic scenario the occurrence of singularity is strictly a four-dimensional 
phenomenon and may be used to hide a nonsmooth transition. 

The second order corrections to the long-wavelength solutions include 
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the effect of spatial gradients. As expected, the picture is symmetric: terms 
which are growing in the pre-big bang phase are decaying in the post-big 
bang phase and vice versa; we might use different terms in different asymp- 
totic regimes to estimate the significance of corrective terms. We conclude 
that in the case of an exotic matter with barotropic index r = or T = 2 
we can achieve the decay of all second order correction terms at the end 
of the pre-big bang stage. However, corrections for the post-big bang stage 
always include terms which are growing in time. Our conclusion indicates 
that a simple model without any dilaton potential or cosmological constant 
leads to a result which can hardly fit our present understandings. Inclusion 
of additional terms in the field equations (e.g. V ((/>), A, or H^ up ) changes 
the form of solutions significally and in that case the analysis of the second 
order is far more complicated. 

ACKNOWLEDGEMENTS 

This work was supported by the Estonian Science Foundation under 
grant No 5026. 

References 

[1] Lifschitz, E. M., and Khalatnikov, I. M., (1963). Adv. Phys. 12, 185. 
[2] Tomita, K., (1993). Phys. Rev. D48, 5634. 

[3] Comer, G. L., Deruelle, N., Langlois, D. and Parry, J. (1994). Phys. 
Rev. D49, 2759. 

[4] Khalatnikov, I. M., Kamenshchik, A. Yu., and Starobinsky, A. A. (2002) 
E-print gr-qc/0204045. 

[5] Tomita, K., and Deruelle, N. (1994). Phys. Rev. D50, 7216. 



17 



[6] Comer, G. L., Deruelle, N., and Langlois, D. (1997). Phys. Rev. D55 
3497. 

[7] Salopek, D. S., and Stewart, J. M. (1992) Class. Quantum Grav. 9, 
1943; Perry, J., Salopek, D. S., and Stewart, J. M. (1994) Phys. Rev. 
D49, 2872. 

[8] Deruelle, N., and Langlois, D. (1995). Phys. Rev. D52, 2007. 

[9] A collection of papers on the pre-big bang scenario is available at home- 
page http://www.to. infn.it/" gasperin/. 

[10] Lidsey, J. E., Wands, D., and Copeland, E. J. (2000). Phys. Rept. 337, 
343. 

[11] Veneziano, G. (1991) Phys. Lett. B265, 287. 

[12] Gasperini, M. and Veneziano, G. (1993) Astropart. Phys. 1, 317. 

[13] Brustein, R. and Veneziano, G. (1994). Phys. Lett. B329, 429. 

[14] Kaloper, M., Madden, R., and Olive, K. A. (1995) Nucl. Phys. B452, 
677. 

[15] Gasperini, M., Ma ggiore, M., and Veneziano, G. (1997). Nucl. Phys. 
B494, 315; Brustein, R., and Madden, R. (1998). Phys. Rev. D57, 712. 

[16] Gasperini, M., Maharana, J., and Veneziano, G. (1996). Nucl. Phys. 
B472, 349; Dabrowski, M. P., and Kiefer, C. (1997). Phys. Lett. B397, 
185. 

[17] Ellis, G. F. R, Roberts, D. C., Solomons, D., and Dunsby, P. K. S. 
(2000) Phys. Rev. D62 084004. 

[18] Infante, M. P. and Sanchez, N. (2000) Current Topics in Astrofun- 
damental Physics: the Cosmic Microwave Background (VII Chalogne 
School, NATO ASI Series C, Kluwer). 



18 



[19] Khoury, J., Ovrut, B. A., Seiberg, N., Steinhardt, P. J., and Turok, N. 
(2002) Phys. Rev. D65, 086007. 

[20] Khoury, J., Ovrut, B. A., Steinhardt, P. J., and Turok, N. (2001) Phys. 
Rev. D64, 123522. 

[21] Steinhardt, P. J. and Turok, N. (2002) Cosmic Evolution in a Cyclic 
Universe, http://feynman.princeton.edu/ steinh . 

[22] Brustein, R., Gasperini,M., Giovannini, M., Mukhanov,V. F., and 
Veneziano, G. (1995) Phys. Rev. D51, 6744. 

[23] Nambu, Y. and Yamaguchi, Y. Y. (1999) Phys. Rev D60, 104011; 
Nambu, Y. (2000) Phys. Rev D62, 104010; Nambu, Y. (2001) Phys. 
Rev D63, 044013. 

[24] Kuusk, P. and Saal, M. (2002) Gen. Rel. Grav. 34, 353. 



19 



